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Abstract 

The  potential  and  current  density  distribution  on  the  electrodes  of  a  lithium-polymer  battery  were  studied  by  using  the  finite  element  method. 
The  effect  of  the  configuration  of  the  electrodes  such  as  the  aspect  ratio  of  the  electrodes  and  the  size  and  placing  of  current  collecting  tabs  as  well 
as  the  discharge  rates  on  the  battery  performance  was  examined  to  enhance  the  uniformity  of  the  utilization  of  the  active  material  of  electrodes. 
The  results  showed  that  the  aspect  ratio  of  the  electrodes  and  the  size  and  placing  of  current  collecting  tabs  have  a  significant  effect  on  the  potential 
and  current  density  distribution  on  the  electrodes  to  influence  the  distribution  of  the  depth  of  discharge  on  the  electrodes,  thus  affecting  the  uniform 
utilization  of  the  active  material  of  electrodes. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

There  is  a  significant  interest  in  the  use  of  batteries  for  hybrid 
electric  vehicle  (HEV)  and  electric  vehicle  (EV).  The  outstand¬ 
ing  characteristics  of  lithium-polymer  batteries  (high  energy 
density,  high  voltage,  low  self-discharge  rate,  and  good  sta¬ 
bility  among  others)  make  them  one  of  the  preferred  choices 
for  such  applications.  However,  much  larger  lithium-polymer 
batteries  than  those  available  in  the  market  for  consumer  elec¬ 
tronics  are  required  for  HEV  and  EV  applications.  The  perfor¬ 
mance  of  a  battery  electrode  is  influenced  by  the  aspect  ratio, 
the  placing  of  current  collecting  tabs,  and  the  total  amount 
of  the  current  flowing  through  an  electrode.  If  an  electrode  is 
not  designed  optimally,  the  potential  and  current  density  will 
be  non-uniformly  distributed,  and  the  utilization  of  the  active 
material  over  the  electrode  will  be  non-uniform.  Accelerated 
degradation  of  the  electrode  may  result  due  to  excessive  local¬ 
ized  utilization  of  the  active  material  on  the  electrode.  That 
effect  becomes  more  pronounced,  as  the  size  of  the  electrode 
becomes  larger.  Therefore,  an  optimum  design  of  the  electrode 
is  pertinent  for  the  production  of  large-scale  lithium-polymer 
batteries. 
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When  scaling  up  a  small-scale  cell  to  a  large-scale  battery, 
mathematical  modeling  plays  an  important  role,  because  nearly 
limitless  design  iterations  can  be  performed  by  using  simula¬ 
tions  [1].  Previous  reviews  of  the  modeling  of  lithium  batteries 
are  given  in  references  [2-5] .  A  one-dimensional  model  assumes 
that  the  gradients  of  the  variables  adopted  in  modeling  are  neg¬ 
ligible  in  the  two  directions  parallel  to  the  current  collectors. 
Such  an  assumption  may  be  valid  for  small-scale  cells.  However, 
that  assumption  may  not  be  justified  for  large-scale  batteries, 
since  the  potential  drop  along  the  current  collector  due  to  ohmic 
drop  may  be  significant  enough  to  affect  the  current  distribution, 
with  a  higher  current  closer  to  the  tabs.  Then,  a  two-  or  three- 
dimensional  model  may  be  desirable  for  large-scale  batteries 
[6-9]. 

In  this  work,  a  two-dimensional  modeling  is  performed  to 
calculate  the  potential  and  current  density  distribution  on  the 
electrodes  of  a  lithium-polymer  battery  comprising  a  LiM^C^ 
cathode,  a  graphite  anode,  and  a  plasticized  electrolyte.  This 
work  adopts  a  relatively  simpler  modeling  approach  by  consid¬ 
ering  only  Ohm’s  law  and  charge  conservation  on  the  electrodes 
based  on  the  simplified  polarization  characteristics  of  the  elec¬ 
trodes  as  compared  to  the  previously  published  papers  by  other 
researchers  [10-15].  The  distribution  of  the  depth  of  discharge 
(DOD)  on  the  electrode  is  predicted  as  a  function  of  discharge 
time  from  the  calculated  potential  and  current  density  distribu¬ 
tion.  Based  on  the  distribution  of  DOD,  the  effects  of  the  aspect 
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ratio  of  the  electrodes,  the  size  and  placing  of  current  collecting 
tabs,  and  discharge  rates  on  the  battery  performance  are  evalu¬ 
ated. 


of  the  positive  and  negative  electrodes,  respectively.  By  substi¬ 
tuting  Eqs.  (3)  and  (4)  into  Eqs.  (1)  and  (2),  the  following  Poisson 
equations  for  Vp  and  Vn  are  obtained: 


2.  Mathematical  model 

A  schematic  diagram  of  the  current  flow  in  the  parallel  plate 
electrodes  of  a  battery  is  shown  in  Fig.  E  The  distance  between 
the  electrodes  is  assumed  to  be  so  small  that  the  current  flow 
between  the  electrodes  is  perpendicular  to  the  electrodes.  From 
the  continuity  of  current  on  the  electrodes,  the  following  equa¬ 
tions  can  be  derived: 

V  -/p  -  J  =  0  in  (1) 

V  •  /n  +  /  =  0  in  (2) 

where  /p  and  in  are  the  linear  current  density  vectors  (current 
per  unit  length  (A  cm-1))  in  the  positive  and  negative  elec¬ 
trodes,  respectively,  and  /  is  the  current  density  (current  per 
unit  area  (A  cm-2))  transferred  through  the  separator  from  the 
negative  electrode  to  the  positive  electrode.  £2p  and  denote 
the  domains  of  the  positive  and  negative  electrodes,  respectively. 


By  Ohm’s  law,  /p  and  in  can  be  written  as 

1 

try  —  ^7  Vp  in  £2p 

(3) 

rP 

1 

i  f)  —  V  E  in 

(4) 

P n 


where  rp  and  rn  are  the  resistances  (£2)  of  the  positive  and  nega¬ 
tive  electrodes,  respectively,  and  Vp  and  Vn  are  the  potentials  (V) 


Current 


Fig.  1.  Schematic  diagram  of  the  current  flow  in  the  parallel  plate  electrodes  of 
a  battery. 


V2VP  =  -rpJ  in  £2p  (5) 

V2Vn  =  +rnJ  in  Qn  (6) 


The  relevant  boundary  conditions  for  Vp  are 
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where  d/dn  denotes  the  gradient  in  the  direction  of  the  outward 
normal  to  the  boundary.  The  first  boundary  condition  (7)  implies 
that  there  is  no  current  flow  through  the  boundary  (rPl )  of  the 
electrode  other  than  the  tab.  The  second  boundary  condition  (8) 
means  that  the  linear  current  density  through  the  tab  (rp2)  of  the 
length  L  (cm)  is  constant  to  the  value  of  Iq/L.  Io  is  the  total  current 
(A)  through  the  tab  in  the  mode  of  constant-current  discharge. 
The  boundary  conditions  for  Vn  are 


—  =  0  on  T 
dn 

yn  =  0  on  rn2 


(9) 

(10) 


The  first  boundary  condition  (9)  implies  the  same  as  in  the  case  of 
Vp.  The  second  boundary  condition  (10)  means  that  the  potential 
at  the  tab  of  the  negative  electrode  is  fixed  to  the  value  of  zero 
as  the  reference  potential. 

The  resistance,  r  (rp  or  rn),  is  calculated  as  follows: 

1 

r  =  -  (11) 

hcSc  +  hQS& 

where  hc  and  he  are  the  thicknesses  (cm)  of  the  current  collector 
and  the  electrode  material,  respectively,  and  Sc  and  SQ  are  the 
electrical  conductivities  (Scm-1)  of  the  current  collector  and 
the  electrode  material,  respectively.  The  parameters  used  in  the 
calculations  of  resistances  for  the  electrodes  are  listed  in  Table  1 . 
The  values  of  electrical  conductivities  of  composite  electrodes 
are  the  same  as  those  used  in  the  references  [14,15].  The  cur¬ 
rent  collectors  of  cathode  and  anode  are  made  of  aluminum  and 
copper,  respectively. 

The  current  density,  /,  of  Eqs.  (5)  and  (6)  is  the  function  of 
the  potential  difference  between  the  positive  and  negative  elec¬ 
trodes,  (Vp  —  Vn).  The  functional  form  depends  on  the  polariza¬ 
tion  characteristics  of  the  electrodes.  In  this  study,  the  following 


Table  1 

Parameters  for  the  electrodes 


Parameter 

LryMn204 

Se  (S  cm- 1 ) 

1.0 

0.038 

h  (|xm) 

85 

140 

Sc  (S  cm- 1 ) 

6.33  x  105 

3.83  x  105 

hc  ([xm) 

10 

20 
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polarization  expression  used  by  Tiedemann  and  Newman  [16] 
and  Newman  and  Tiedemann  [9]  was  adopted 


J  =  Y(V p  -Vn~U)  (12) 

where  Y  and  U  are  the  fitting  parameters.  As  suggested  by  Gu 
[17],  U  and  Y  were  expressed  as  the  following  functions  of  the 
depth  of  discharge: 

U  =  ao  +  <21  (DOD)  +  a2(DOD)2  +  <23  (DOD)3  (13) 

Y  =  <24  +  <25  (DOD)  +  <36(DOD)2  (14) 


electrode  can  be  calculated  from  the  distribution  of  J  as 


DOD  = 


(15) 


where  tis  the  discharge  time  (s)  and  Qj  is  the  theoretical  capacity 
per  unit  area  (Ah  cm-2)  of  the  electrodes.  The  uniformity  index, 
UI,  of  DOD  is  defined  as 


UI  = 


DODmax  DODmin 

2DODavg 


(16) 


where  a^-a 6  are  the  constants  to  be  determined  by  experiments. 

By  solving  the  equations  listed  previously,  the  distribution 
of  the  current  density,  /,  on  the  electrodes  can  be  obtained  as  a 
function  of  the  position  on  the  electrode  and  the  time.  Therefore, 
DOD  varies  along  with  the  position  on  the  electrode  and  the 
time  elapsed  during  discharge.  The  distribution  of  DOD  on  the 


where  DODmax,  DODmin,  and  DODavg  are  the  maximum,  mini¬ 
mum,  and  average  values  of  DOD  on  the  electrodes,  respectively. 
The  lower  the  value  of  UI,  the  better  the  uniformity  of  DOD. 
Hence,  a  way  to  reduce  UI  should  be  devised  in  order  to  improve 
the  uniformity  of  the  utilization  of  the  active  material  of  the  elec¬ 
trodes  and  to  lengthen  the  life  cycle  of  the  battery. 
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Fig.  2.  Schematic  diagrams  of  the  electrode  shapes  for:  (a)  types  A-F,  (b)  type  G,  and  (c)  type  H. 
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Capacity  (mAh) 

Fig.  3.  Comparison  between  experimental  and  modeling  discharge  curves  at 
discharge  rates  of  1C,  2 C,  3C,  5C,  8 C,  and  IOC.  Solid  lines  are  experimental 
data  and  lines  with  squares  are  modeling  results  based  on  the  Unite  element 
method. 


3.  Results  and  discussion 

The  solutions  to  the  governing  Eqs.  (5)  and  (6)  subject  to 
the  associated  boundary  conditions  (7)— (10)  were  obtained  by 
using  the  finite  element  method.  Numerical  simulations  were 
performed  for  the  electrodes  of  eight  different  shapes  of  a 
lithium-polymer  battery  having  a  nominal  capacity  of  5  Ah. 
The  electrodes  of  eight  different  shapes  are  referred  to  as  the 
electrodes  of  types  A-H.  The  schematic  diagrams  of  the  elec¬ 
trode  shapes  for  types  A-F,  type  G,  and  type  H  are  illustrated 
in  Fig.  2(a-c),  respectively.  The  dimensions  of  the  electrodes 
and  the  size  and  placing  of  current  collecting  tabs  are  listed 
in  Table  2.  In  order  to  test  the  validity  of  the  modeling,  the 
calculated  discharge  curves  based  on  modeling  are  compared 
with  the  experimental  data  in  Fig.  3.  The  experiments  were  per¬ 
formed  at  room  temperature  by  using  the  5  Ah  battery  fabricated 
by  VK  Corporation  with  the  electrodes  of  type  A,  of  which 
the  dimensions  of  the  electrodes  and  the  positions  of  the  tabs 


VP 

4.08698 

4.08694 

4.08687 

4.08680 

4.08672 

4.08665 

4.08658 

4.08651 

4.08644 

4.08640 


VP 

3.86951 

3.86948 

3.86941 

3.86933 

3.86926 

3.86919 

3.86912 

3.86904 

3.86897 

3.86894 


VP 

2.96772 

2.96769 

2.96761 

2.96754 

2.96747 

2.96740 
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2.96725 

2.96718 

2.96714 


8 
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Fig.  4.  Change  of  the  distributions  of  potential  on  the  positive  electrode  of  type  A  during  the  discharge  with  1C  rate.  Discharge  times  are:  (a)  1  min,  (b)  30  min,  and 
(c)  59  min. 
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Table  2 

Dimensions  of  the  electrodes  and  the  size  and  placing  of  current  collecting  tabs 
(unit:  mm) 
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W 

H 

D 

L 

Type  A 

86 

149 

3 

35 

Type  B 

149 

86 

3 

35 

Type  C 

113.2 

113.2 

3 

35 

Type  D 

86 

149 

3 

18 

Type  E 

86 

149 

25 

35 

Type  F 

86 

149 

34 

18 

Type  G 

Refer  to  Fig.  2(b) 

Type  H 

Refer  to  Fig.  2(c) 

are  shown  in  Table  2  and  Fig.  2(a).  At  various  discharge  rates 
from  1C  to  IOC,  the  experimental  discharge  curves  are  in  good 
agreement  with  the  modeling  results  based  on  the  finite  element 
method. 

In  Fig.  4,  the  change  of  the  distributions  of  potential  on  the 
positive  electrode  of  type  A  during  the  discharge  with  1 C  rate 
is  shown.  The  potential  gradient  is  seen  to  be  most  severe  in 
the  region  where  the  tab  is  attached  to  the  current  collector. 
This  is  because  all  the  current  flows  through  the  conducting 
current  collector  into  the  tab  from  the  entire  electrode  plate. 
Note  that  the  values  of  the  potential  of  Fig.  4(a-c)  are  around 
4.09,  3.87,  and  2.97  V,  respectively.  Although  the  appearances 
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-0.00059 

-0.00063 


Fig.  5.  Distribution  of  potential  on  the  negative  electrode  of  type  A  for  the 
discharge  rate  of  1 C  at  the  discharge  time  of  1  min. 


of  the  potential  distribution  are  similar,  the  values  of  the  poten¬ 
tial  at  the  beginning,  in  the  middle,  and  at  the  end  of  discharge 
decrease  complying  with  the  discharge  characteristics  of  the  bat¬ 
tery.  In  Fig.  5,  only  the  distribution  of  potential  on  the  negative 
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Fig.  6.  Change  of  the  distributions  of  current  density  transferred  through  the  separator  from  the  negative  electrode  to  the  positive  electrode  of  type  A  during  the 
discharge  with  1C  rate.  Discharge  times  are:  (a)  1  min,  (b)  30  min,  and  (c)  59  min. 
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Fig.  7.  UI  of  DOD  of  the  electrode  of  type  A  as  a  function  of  the  average  DOD 
for  the  discharge  rates  of  1C,  2C,  3C,  5C,  and  IOC. 

electrode  of  type  A  at  the  discharge  time  of  1  min  is  shown, 
because  the  values  of  the  potential  and  the  appearances  of  the 
potential  distribution  on  the  negative  electrode  do  not  change 
according  to  the  state  of  discharge.  Again,  the  potential  gradi¬ 
ent  is  highest  at  the  region  near  tab,  because  all  the  current  has 
to  flow  from  the  tab  through  the  entire  electrode  plate.  Fig.  6 
shows  the  change  of  the  non-uniform  distributions  of  current 
density  transferred  from  the  negative  electrode  to  the  positive 
electrode  of  type  A  during  the  discharge  with  1C  rate.  At  the 
initial  stage  of  discharge  (Fig.  6(a)),  current  density  is  higher 
in  the  region  near  the  tabs  than  in  the  region  distant  from  the 
tabs.  However,  after  the  midway  stage  of  discharge  (Fig.  6(b 
and  c)),  current  density  is  lower  in  the  region  near  the  tabs  than 
in  the  region  distant  from  the  tabs.  This  is  because  the  less  uti¬ 
lized  region  distant  from  the  tabs  becomes  more  favorable  for 
reaction  after  the  midway  stage  of  discharge  than  the  region 
near  the  tabs  which  is  heavily  utilized  at  the  initial  stage  of  dis¬ 
charge.  A  ‘deformation’  in  the  current  density  distribution  at  the 
upper  center  of  the  electrode  observed  in  Fig.  6  may  be  due  to 
the  abrupt  change  of  the  directions  of  the  current  flow  near  the 
corners  of  the  tabs  of  the  positive  and  negative  electrodes.  The 
distributions  of  DOD  on  the  electrodes  can  be  obtained  from 
the  integration  of  the  current  density  with  time  shown  in  Fig.  6. 
Although  the  DOD  is  higher  in  the  region  near  the  tabs  than  in 
the  region  distant  from  the  tabs  throughout  the  whole  stages  of 
discharge,  the  difference  between  the  maximum  and  minimum 
values  of  DOD  decreases  as  the  discharge  proceeds,  since  the 
reaction  rate  tends  to  be  higher  in  the  region  where  the  active 
material  is  less  utilized.  This  is  reflected  in  the  change  of  UI 
of  DOD  with  time.  At  the  initial  stage  of  discharge,  the  value 
of  UI  of  DOD  is  high,  but  it  decreases  as  the  discharge  goes 
on.  In  Fig.  7,  the  UI  of  DOD  for  the  discharge  rates  of  1C,  2C, 
3 C,  5C,  and  IOC  are  shown  as  a  function  of  the  average  DOD 
over  the  entire  electrode  area.  At  all  discharge  rates,  the  UI  of 
DOD  decreases  as  average  DOD  increases.  The  UI  of  DOD 
maintains  a  higher  value  for  a  higher  discharge  rate  than  for  a 
lower  discharge  rate  throughout  the  whole  stages  of  discharge. 
This  trend  illustrates  that  the  non-uniformity  of  the  utilization  of 
the  active  material  of  electrodes  becomes  more  pronounced  at 
high  discharge  rates.  Therefore,  when  the  high  discharge  rates 


Fig.  8.  UFs  of  DOD  of  the  electrodes  of  types  A-C  as  a  function  of  the  average 
DOD  for  the  discharge  rate  of  1C. 

are  repeated,  there  is  a  high  chance  that  only  a  small  fraction 
of  electrode  will  be  utilized  excessively,  thus  reducing  the  life 
cycle  of  the  battery. 

To  investigate  the  effect  of  the  aspect  ratio  of  the  electrodes 
on  the  uniformity  of  the  utilization  of  the  active  material  of 
electrodes,  the  UI’s  of  DOD  for  the  electrodes  of  types  A-C  are 
shown  as  a  function  of  the  average  DOD  for  the  discharge  rate 
of  1C  in  Fig.  8.  It  is  assumed  for  the  sake  of  modeling  that  the 
electrodes  A-C  have  the  same  surface.  In  terms  of  the  UI  of 
DOD,  the  electrode  of  type  B  ( W>H)  has  advantage  over  the 
electrode  of  type  A  (W<H)  or  the  electrode  of  type  C  (W=H), 
where  W  and  H  denote  the  width  and  height  of  the  electrode, 
respectively.  In  Fig.  9,  the  UI’s  of  DOD  of  the  electrodes  of 
types  A,  D-G,  and  H  are  plotted  as  a  function  of  the  average 
DOD  for  the  discharge  rate  of  1 C  in  order  to  study  the  effect  of 
the  size  and  placing  of  current  collecting  tabs  on  the  uniformity 
of  the  utilization  of  the  active  material  of  electrodes.  This  plot 
clearly  shows  which  type  of  electrode  is  favorable  in  terms  of 
the  UI  of  DOD.  From  the  plots  of  the  UI  of  DOD  of  the  types 
A  and  D,  it  can  be  deduced  that  the  enlargement  of  the  tab  size 


Fig.  9.  UFs  of  DOD  of  the  electrodes  of  types  A,  D-G,  and  H  as  a  function  of 
the  average  DOD  for  the  discharge  rate  of  1 C. 
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is  favorable  to  the  uniform  utilization  of  active  material.  Based 
on  the  comparisons  between  the  types  A  and  E  and  between  the 
types  D  and  F,  the  reduction  of  the  distance  between  the  tabs  is 
unfavorable  for  the  cases  investigated.  As  a  major  modification 
of  the  placing  of  tabs,  the  UI’s  of  DOD  as  a  function  of  the 
average  DOD  are  calculated  for  the  types  H  and  G,  where  the 
tabs  are  on  the  opposing  sides  of  a  rectangular  electrode.  The 
electrodes  of  types  H  and  G  is  more  favorable  than  the  others  in 
terms  of  the  UI  of  DOD.  By  the  comparison  between  the  types 
H  and  G,  it  can  be  noticed  that  placing  the  tabs  diagonally  is 
favorable,  although  the  difference  between  the  UI’s  of  DOD  of 
types  H  and  G  is  minute. 

4.  Conclusions 

A  mathematical  procedure  was  developed  to  study  the  poten¬ 
tial  and  current  density  distribution  on  the  electrodes  for  the 
scale-up  of  the  electrode  of  a  lithium-polymer  battery.  The  dis¬ 
tribution  of  the  depth  of  discharge  on  the  electrode  was  predicted 
as  a  function  of  discharge  time  from  the  calculated  potential  and 
current  density  distribution  by  using  the  finite  element  method. 
A  quantitative  examination  was  made  of  the  effect  of  varying 
the  configuration  of  the  electrodes  such  as  the  aspect  ratio  of 
the  electrodes  and  the  size  and  placing  of  current  collecting  tabs 
as  well  as  the  discharge  rates  on  the  battery  performance,  in 
order  to  enhance  the  uniformity  of  the  utilization  of  the  active 
material  of  electrodes.  The  results  showed  that  the  aspect  ratio 
of  the  electrodes  and  the  size  and  placing  of  current  collecting 
tabs  have  a  significant  effect  on  the  potential  and  current  density 
distribution  on  the  electrodes  to  influence  the  DOD  distribution 
on  the  electrodes,  thus  affecting  the  uniform  utilization  of  the 
active  material  of  electrodes.  The  methodology  presented  in  this 
study,  used  to  evaluate  quantitatively  the  uniformity  of  the  uti¬ 
lization  of  the  active  material  of  electrodes,  should  contribute  to 
modifying  the  design  of  the  electrode  configuration. 
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